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1.  Introduction 


Two-dimensional  (2-D)  operations,  such  as  rotations,  translations,  and  intersections,  are  tools 
that  are  essential  for  many  types  of  scientific  modeling.  However,  the  C++  programming 
language  does  not  natively  perform  them.  This  report  documents  a  set  of  functions,  written  in 
C++,  that  can  be  used  to  perform  2-D  rotations,  translations,  and  intersections.  All  of  the 
functions  have  been  grouped  into  the  y2DOps  namespace,  which  is  summarized  at  the  end  of  this 
report. 

The  functions  that  are  presented  in  this  report  are  special  cases  of  more  general  three- 
dimensional  (3-D)  functions.1  Compared  to  the  3-D  functions,  the  2-D  functions  provide  simpler 
interfaces  and  faster  calculations. 


2.  Translation  of  a  Point  in  Space 


2.1  Derivation 

Let  the  position  vector  p  represent  an  arbitrary  point  in  a  plane,  where 

p  =  pxx+  pj  .  (1) 

Furthermore,  let  d  represent  a  displacement  vector,  where 

d  =  dxx  +  dyy .  (2) 

If  p  is  used  to  represent  p  after  it  has  been  translated,  then 

P  =  (p,  +  d,  )-x  +  (py  +  dy  )>’  •  (3) 


2.2  C++  Implementation 
Translate2D()  Code 


inline  void  Translate2D(//<==========================PERF0RMS  A  2D  TRANSLATION 

double  p[2],//< . COORDINATES  TO  TRANSLATE  (MODIFIED  BY  THIS  FUNCTION) 

const  double  d  [2]  ){//< . DISPLACEMENT  VECTOR 

p[0]+=d[0]  ,  p[l]+=d[l]; 

^  j  j  ivm/v/vY  AG  E  N  AU  T  ^GMA  XL*  C  LAS  T  ^  LJ  P  DAT  E  0  2  MAY  201 


Aager,  R.  J.  Three-Dimensional  Translations,  Rotations,  and  Intersections  Using  C+  +  ;  U.S.  Army  Research  Laboratory: 
Aberdeen  Proving  Ground,  MD,  2013,  in  press. 
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Translate2D()  Parameters 

p  p  is  a  two-element  array  that  stores  the  position  vector  that  is  described  by  equation  1 
(P={PV’  Py }  )•  Note  that  p  is  modified  by  the  Translate()  function,  as  described  by 

equation  3. 

d  d  is  a  two-element  array  that  stores  the  displacement  vector  that  is  described  by 
equation  2  (d=  {dx,d  }).  d  determines  the  amount  and  direction  by  which  p  is 

translated. 

Translate2D()  Example 

Figure  1  shows  point  p  being  translated  to  a  new  position  (  p'  )  by  displacement  vector  d . 


Figure  1 .  Translate()  example. 

Let  p  =  {3,l}  and  d  =  {-  2,l}.  Point  p'  can  be  found  by  using  the  Translate2D()  function,  as 
shown  in  the  following  sample  code. 


#include  <cstdio>// . printfQ 

#include  "y_2d_ops .  h"// . Translate2D() 

int  main(){ 

double  p[2]={3,l}; 
double  d[2]={-2,l}; 


y2D0ps: :Translate2D(p,d); 
printf("p[0]=%f,  p[l]=%f\n"jP[0],p[l])j 
^  j  I  <v»vivfv/vivY  AG  E  N  AUT  (SGMA  XL*  C  OFlrwrw,'jrurw,v,v,>jrv,v,v,v,>j,v/v,v,v,v,vrv,v,v,v,,v,Nj  LAS  T^U  P  DAT  2F1AY  201. 


OUTPUT: 

p [0] =1 . 000000J  p [ 1 ] =2 . 000000 
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3.  Rotation  of  a  Point  About  an  Arbitrarily  Positioned  Axis 


3.1  Derivation 

Suppose  that  the  unit  vector  v  is  used  to  define  an  arbitrary  axis  about  which  a  point  in  space 
will  be  rotated. 

v  =  vxx  +  vyy  +  vzz .  (4) 

Rodrigues’s  rotation  formula  can  be  used  to  construct  a  rotation  matrix,2  R  ,  that  can  be  used  to 
perform  a  rotation  about  v  by  an  angle  0 .  The  direction  of  the  rotation  can  be  determined  by 
using  the  right-hand-thumb  rule  (when  the  right  thumb  is  pointed  in  the  direction  of  v ,  the  curled 
fingers  of  the  right  hand  will  point  in  the  direction  of  the  rotation). 


where 


~Cg)+Cg 

vxvy(l~ca)-vzse 

vxvf~ce)+  vyse 

-ce)+vzss 

v2y$-ce)+ce 

(5) 

-Co)~VyS0 

Vz(1_C*)+VA 

V2z(l~Cg)+Cg 

cg  =  cos(0) 

(6) 

Sg  =  sin(6>) . 

(7) 

For  the  2-D  case,  assume  that  v  points  in  the  positive  z  direction.  Then 

v  =  0 ,  v  =  0 ,  and  v  =  1 . 


This  greatly  simplifies  equation  5: 


ce  -sg  0 


R  =  |  Sg  Cg  0  • 

0  0  1 


Substituting  equations  6  and  7  into  equation  5,  then  converting  to  2D, 

cos  6  -  sin  0 
R  = 

sin  6  cos  6 


zMason,  M.  T.  Mechanics  of  Robotic  Manipulation ;  Massachusetts  Institute  of  Technology  Press:  Cambridge,  MA,  2001, 
(p  46,  equation  3.26). 
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Let  the  position  vector  p  locate  an  arbitrary  point  in  a  plane,  where 

p  =  pxx  +  pyy  .  (11) 

Let  the  position  vector  3  locate  the  origin  of  the  rotation  axis  defined  by  v . 

d  =  oxx  +  oyy.  (12) 

The  translation-rotation-translation  sequence  described  by  equation  13  can  be  used  to  find  p' , 
where  p'  is  used  to  represent  p  after  it  has  been  rotated  about  v . 

p'  =  R(p-o)+o.  (13) 


3.2  C++  Implementation 

Two  functions  are  used  to  perform  2-D  rotations.  The  first  function,  RMatrix2D(),  calculates  the 
rotation  matrix  that  is  presented  in  equation  10.  The  second  function,  Rotate2D(),  performs  the 
rotation  that  is  presented  in  equation  13.  Breaking  the  calculation  into  two  functions  allows 
functions  that  rotate  objects  containing  more  than  one  point  to  be  written  in  a  manner  that 
doesn’t  sacrifice  performance. 

Rmatrix2D()  Code 

inline  void  RMatrix2D(//<======================CALCULATES  A  2D  ROTATION  MATRIX 

double  R[4] ,  / / < . ROTATION  MATRIX  (CALCULATED  BY  THIS  FUNCTION) 

double  rads){//< . -THE  ANGLE  OF  THE  ROTATION  (CCW  IS  POSITIVE) 

R[0]=R[3]=cos(rads)  ,  R[l]=-(R[2]=sin(rads)); 

}//  rsjrsjrsjtsj  YAGENAUT@GMAIL.COM - LAST-UPDATED-02MAY2013 - 


Rmatrix2D()  Parameters 

R  R  is  a  four-element  array  that  stores  the  rotation  matrix  that  is  described  by  equation  10 
(R=  { R() 0,R01,R10,R11}).  Note  that  R  is  modified  by  the  Rmatrix2D()  function.  R  is 

intended  to  be  used  as  the  third  argument  of  the  Rotate2D()  function. 

rads  rads  is  used  to  represent  the  angle  (in  radians)  of  the  rotation.  The  direction  of  the 
rotation  is  counterclockwise  (see  figure  2). 

Rotate2D()  Code 

inline  void  Rotate2D(//<===================================PERF0RMS  A  ROTATION 

double  p [2 ],//< . COORDINATES  TO  ROTATE  (MODIFIED  BY  THIS  FUNCTION) 

const  double  o[2]J//<-- . -THE  ORIGIN  OF  THE  AXIS  OF  ROTATION 

const  double  R[4]){//< . A  ROTATION  MATRIX  (FROM  RMatrix2D()) 

double  t0=p[0]-o[0]  ,  tl=p[l] -o[l] ; 

p[0]=R[0]*t0+R[l]*tl+o[0]  ,  P [ 1 ] =R [ 2 ] *t0+R [ 3 ] *tl+o [ 1 ] ; 

}//  r\j  r\j  r\j  r\j  YAGENAUT (5)GMAI  L  .  L AST ~U P DAT E D~0 2MAY 201 
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Rotate2D()  Parameters 

p  p  is  a  two-element  array  that  stores  the  position  vector  that  is  described  by  equation  1 1 
(P={PV’  Py }  )•  Note  that  p  is  modified  by  the  Rotate2D()  function,  as  described  by 

equation  13. 

o  o  is  a  two-element  array  that  stores  the  position  vector  that  is  described  by  equation  12 
(o=  {o x ,  oy } ).  o  is  the  point  about  which  p  is  rotated. 

R  R  is  a  rotation  matrix  that  has  been  precalculated  using  the  RMatrix2D()  function. 

Rotate2D()  Example 

Figure  2  shows  point  p  being  rotated  about  o  to  a  new  position  ( p' ). 


Figure  2.  Rotate()  example. 

Let  p  =  {3,l}  and  o  =  {2,l}.  Furthermore,  let  the  angle  of  rotation  be  n  12.  Point  p  can  be  found 
by  first  using  the  RMatrix2D()  function  to  calculate  a  rotation  matrix,  then  using  the  Rotate2D() 
function  to  perform  the  rotation. 


#include  <cstdio>// . printfQ 

#include  "y_2d_ops.h 
int  main(){ 

double  p[2]={3,l}; 
double  o[2]={2,l}; 

double  R[4] ;  /*<-*/y2D0ps : : RMatrix2D(R, 3. 14159265358979/2) ; II...  rotation  matrix 

y2D0ps: : Rotate2D(p, o, R) ; 

printf ("p[0]=%f ,  p[l]=%f\n"Jp[0],p[l])J 

J I  j  »v»MivfN;/vivY  AG  E  N  AUT  ^)GMA  XL*  C  LAST  ^  LJ  P  DAT  E  D  ^  0  2MAY  201.3  fv/viviv/N/iv 


OUTPUT: 

p [0] =2 . 000000,  p [ 1 ] =2 . 000000 
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4.  Intersection  Between  Two  Lines 


4.1  Derivation 

Suppose  that  line  A  passes  through  the  points  A()  and  A,  where 

A)  =  A), A + A.v.'v  and  A  =  A,A + A.A  •  (14) 

Let  pA  represent  a  point  that  lies  on  A  .  A0  and  A,  can  be  used  to  construct  a  parametric 
equation  for  pA  as  a  function  of  the  parameter  t0 : 

Pa  =  A  +  (A  ~A>)to- 

The  parameter  t0  represents  the  scaled  distance  from  A0  to  A  along  A  .  Thus,  if  t()  =  0 ,  pA  is 
located  at  A0 .  If  t0  =  1 ,  pA  is  located  at  A, . 

Similarly,  suppose  that  line  B  passes  through  the  points  B0  and  B{  where 

B0  =  B0xx  +  B0y y  and  Bx  =  Blxx  +  Blyy.  (16) 

Let  pB  represent  a  point  that  lies  on  B  .  B()  and  5,  can  be  used  to  construct  a  parametric 
equation  for  pB  as  a  function  of  the  parameter  t{ : 

pB=B0  +  (Bl-B0)t1.  (17) 

Figure  3  presents  an  image  of  lines  A  and  B  for  the  case  where  they  intersect. 
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The  point  of  intersection  between  A  and  B  occurs  where  pA  is  equal  to  pB .  Thus, 

A)  +  (a  ~  A))fo  =  b0  +  {b1  —  B0^t1. 

Rearranging  terms, 

A, -Bo  =  (A)  —  A)^o  +  (®i  ~ 


This  can  be  written  in  matrix  form  as 


Au  B0  x 

1 

1 

Au  —  B0  y  - 

1 

1 

Au 

A),y 


B\,x  B0x  t0 

B\y  —  B0y  tl 


Solving  for  t  , 


^0 

i 

cA 

i 

P 

u 

Bu 

-V 

-l 

i 

.> 

i 

Co 

© 

_ 1 

_A  ,y  A,y 

K 

-V 

A,y  B()  y  _ 

(18) 


(19) 


(20) 


(21) 


Recall  that  the  vector  t  contains  the  parameters  from  equations  15  and  17.  Thus,  once  t  is 
known,  t0  can  be  substituted  into  equation  15  to  find  the  point  of  intersection  between  A  and  B  . 

Note  that  if  the  two-by-two  matrix  defined  in  equation  21  is  noninvertible,  then  A  is  parallel  to 
B . 
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Because  t0  and  tl  are  defined  to  be  scaled  distances  from  A(]  to  Al  and  B0  to  B, ,  respectively, 
if  the  conditions  presented  in  equation  22  are  met,  then  the  point  of  intersection  lies  between  A0 
and  Aj  and  between  B0  and  B, ,  as  shown  in  figure  3. 

0  <  tx  <1  and  0  <  t2  <  1 .  (22) 


4.2  C++  Implementation 

Two  functions  are  used  to  find  line-line  intersections.  The  first  function,  IParameters2D(), 
calculates  a  two-element  array  that  is  the  solution  to  equation  21.  The  second  function, 
Intersect2D(),  calculates  the  point  of  intersection  between  two  lines. 

Because  there  is  a  chance  that  the  two-by-two  matrix  shown  in  equation  21  will  be  singular,  a 
Boolean  that  indicates  whether  or  not  a  solution  is  valid  is  returned  by  the  IParameters2D() 
function. 

IParameters2D()  Code 


inline  bool  IParameters2D(//<============PARAMETERS  FOR  LINE-LINE  INTERSECTION 

double  t[2],//< . INTERSECTION  PARAMETERS  (CALCULATED  BY  THIS  FUNCTION) 

const  double  A[4],//< . - . -LINE  A  {A0X, A0Y, A1X,A1Y} 

const  double  B[4],//< . -LINE  B  {B0XJB)YJB1XJB1Y} 

double  e=lE-9){//< - CUTOFF  VALUE  FOR  DETERMINING  IF  A  AND  B  ARE  PARALLEL 

double  a=A[0]-A[2]  ,  b=B[2] -B[0] , // . 2x2  matrix 

c=A[l]-A[3]  ,  d=B[3] - B [ 1 ] ; //  elements 

double  D=a*d-b*c;// . determinant  of  2x2  matrix 

if (fabs(D)<e)return  false;//. .. .=>  A  &  B  are  parallel  (and  t  is  meaningless) 
double  f=A[0] -B[0]  ,  g=A[l]-B[l]; 
t[0]=(d*f -b*g)/D  ,  t[l]=(-c*f+a*g)/D; 

return  true;// . =  >  A  &  B  are  not  parallel 

^  j  j  r\ir\ir\>r\>  Y  /\(j  ^  N  AU  T  (fl)GMA  XL*  C  ON- ^  LAST  ^  LJ  P  DAT  E  D  ^  0  2  MAY  2013 


IParameters2D()  Parameters 

R  R  is  a  rotation  matrix  that  has  been  precalculated  using  the  RMatrix2D()  function. 

t  t  is  a  two-element  array  that  stores  the  parameters  described  in  equation  21  (t=  [t0 ,/,}). 
Note  that  t  is  modified  by  the  IParameters2D()  function. 

A  A  is  a  four-element  array  that  stores  the  line  that  is  defined  by  equation  14 

(A=  { A),*  ’  A),y  ’  A,JC  ’  )  )• 

B  B  is  a  four-element  array  that  stores  the  line  that  is  defined  by  equation  16 
(B=  {B0  x ,  B0  y ,  BX  x ,  BX  y } ). 

e  e  is  the  cutoff  value  for  testing  whether  or  not  A  and  B  are  parallel.  If  the  determinant 
of  the  matrix  in  equation  21  is  less  than  e,  then  A  and  B  are  considered  to  be  parallel. 
The  default  value  of  e  is  10  9 . 
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IParameters2D()  Return  Value 

IParameters2D()  returns  false  if  A  is  parallel  to  B.  A  return  value  of  false  indicates  that  t  has  not 
been  calculated  and,  thus,  should  not  be  passed  to  the  Intersect2D()  function. 

Intersect2D()  Code 


inline  bool  Intersect2D(//<============CALCULATES  LINE-LINE  INTERSECTION  POINT 

double  x[2 ],//< . POINT  OF  INTERSECTION  (CALCULATED  BY  THIS  FUNCITON) 

const  double  t [2] , //< . INTERSECTION  PARAMETERS  (FROM  IParameters2D( ) ) 

const  double  A[4]  ){//< . -LINE  A  {A0X, A0Y, A1X.A1Y} 

x[0]=A[0]+t[0]*(A[2]-A[0])  ,  x[l]=A[l]+t[0]*(A[3]-A[l])j 
return  t [0] >0&&t[0]<l&&t[l] >0&&t[l]<l; 

}//  YAGENAUT@GMAIL.COM - LAST~UPDATED~02MAY2013 - 


Intersect2D()  Parameters 

x  x  is  a  two-element  array  that  stores  the  point  of  intersection  between  lines  A  and  B. 
Note  that  x  is  modified  by  the  Intersect2D()  function. 

t  t  is  a  parameter  list  that  has  been  precalculated  using  the  IParameters2D()  function. 

A  A  is  a  two-element  by  two-element  array  that  stores  the  line  that  is  defined  by  equation 
14  (A=  [A0x ,  A0  y ,  Al  x ,  Al  y } ). 

Intersect2D()  Return  Value 

Intersect2D()  returns  true  if  line  segment  A  intersects  line  segment  B. 

Intersect2D()  Example 

Figure  4  shows  intersecting  line-segments  A  and  B  . 


Figure  4.  Intersect2D()  example. 
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Let  j4q  =  {2,0},  A,  =  {2,2},  B0  =  {l,l},  and  Bx  =  {3,l}.  The  point  of  intersection  can  be  found  by 
first  calling  the  IParameters2D()  function,  then  using  the  result  in  the  Intersect2D()  function. 

#include  <cstdio>// . printfQ 

#include  "y_2d_ops . h" 
int  main(){ 

double  A[4]={2,0  ,  2,2};// . a  line  segment 

double  B[4]={1,1  ,  3,1};// . a  line  segment 

double  t [2] ; /*<-*/y2D0ps : : IParameters2D(t , A, B) ; // . intersection  parameters 

double  p[2] ; /*<-*/y2D0ps : : Intersect2D(p, t ,A); // . point  of  intersection 

printf("p[0]=%f,  p[l]=%f \n", p[0] ,p[l] ) ; 

^  j  j  m/vivm/v/vY  AG  E  N  AUT  (SGMAI  | #  G  LAS  T/v,U  P  DAT  E  D  ^  0  2NAY  201  3'n"n"'"v'v"n' 


OUTPUT: 


p [0] =2 . 000000,  p[l]=l. 000000 


5.  Summary 


A  summary  sheet  is  provided  at  the  end  of  this  report.  It  presents  the  y2DOps  namespace,  which 
contains  the  five  functions  that  are  described  in  detail  in  sections  2,  3,  and  4.  Also  presented  are 
two  examples  that  demonstrate  the  versatility  of  the  functions  described  in  this  report.  The  first 
uses  the  Rotate2D()  function  to  calculate  a  set  of  points  that  defines  a  simple  orbit  of  a  moon 
around  a  planet,  which  in  turn  is  in  orbit  around  a  star.  The  second  uses  the  Intersect2D() 
function  to  draw  a  four- sided  spiral.  Both  functions  create  text  files  that  contain  all  of  the 
information  needed  to  create  the  two  images  presented  in  the  summary  sheet. 
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Y2DOps  Summary 


y_2d_ops . h 


#ifndef  Y_2D_0PS_H. 

#define  Y_2D_0PS_H. 

#include  <cmath>// 
namespace  y2D0ps{// 

inline  void  Translate2D(//<; 

double  p[2],//< . 

const  double  d[2]){//< 
p[0]+=d[0]  ,  p[l]+=d[l] ; 

}// - YAGENAUT@GMAIL.  COM - LAST~UPDATED~02MAY2013 - 

inline  void  §0tate2jl(//<===================================PERFORMS  A  ROTATION 

double  p[2]j//< - COORDINATES  TO  ROTATE  (MODIFIED  BY  THIS  FUNCTION) 

const  double  o[2],//< - - - . THE  ORIGIN  OF  THE  AXIS  OF  ROTATION 

const  double  R[4]){//< . . A  ROTATION  MATRIX  (FROM  RMatrix2D()) 

double  t0=p[0]-o[0]  ,  tl=p[l] -o[l] ; 

p[0]=R[0]*t0+R[l]*tl+o[0]  ,  p[l]=R[2] *t0+R[3]*tl+o[l] ; 

}// - YAGENAUT@GMAIL.  COM - LAST~UPDATED~02MAY2013 - 

inline  void  RMatrix2D(//<======================CALCULATES  A  2D  ROTATION  MATRIX 

double  R[4]j//< . --ROTATION  MATRIX  (CALCULATED  BY  THIS  FUNCTION) 

double  rads){//< . . THE  ANGLE  OF  THE  ROTATION  (CCW  IS  POSITIVE) 

R[0]=R[3]=cos(nads)  ,  R[l]=-(R[2]=sin(nads) ) ; 

}// - YAGENAUT@GMAIL.COM - LAST~UPDATED~02MAY2013 - 

inline  bool  Intersect2D(//<============CALCULATES  LINE-LINE  INTERSECTION  POINT 

double  x[2] ,//<--- . POINT  OF  INTERSECTION  (CALCULATED  BY  THIS  FUNCITON) 

const  double  t[2],//< - INTERSECTION  PARAMETERS  (FROM  IParameters2D( ) ) 

const  double  A[4]){//< - LINE  A  {  A0X, A0Y, A1X, A1Y} 

x[0]=A[0]+t[0]*(A[2]-A[0])  ,  x[l]=A[l]+t[0]*(A[3]-A[l])j 
return  t [0] >0&&t [0] <l&&t [1] >0&&t[l]<l; 

}// - YAGENAUT@GMAIL.COM - LAST~UPDATED~02MAY2013 - 

inline  bool  IParameter.s2§(//<============PARAMETERS  FOR  LINE -LINE  INTERSECTION 

double  t[2],//< . INTERSECTION  PARAMETERS  (CALCULATED  BY  THIS  FUNCTION) 

const  double  A[4],//< - - - LINE  A  { A0X , A0Y , A1X , A1 Y } 

const  double  B[4],//< - LINE  B  {B0X, B)Y, BIX, B1Y} 

double  e=lE-9){//< - CUTOFF  VALUE  FOR  DETERMINING  IF  A  AND  B  ARE  PARALLEL 

double  a=A[0]-A[2]  ,  b=B[2] -B[0],// . 2x2  matrix 

c=A[l]-A[3]  ,  d=B[3] -B[l];//  elements 

double  D=a*d-b*c;// . determinant  of  2x2  matrix 

if (fabs(D)<e)return  false;//. .. .=>  A  &  B  are  parallel  (and  t  is  meaningless) 
double  f=A[0]-B[0]  ,  g=A[l]-B[l]; 
t[0]=(d*f-b*g)/D  ,  t[l]=(-c*f+a*g)/D; 

return  true;// . =  >  A  &  B  are  not  parallel 

}// - YAGENAUT@GMAIL .  COM - LAST~UPDATED~02MAY2013 - 

} / / @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ 

#endif 


EXAMPLE  1 


#include  <cstdio>// . FILE, f reopen ( ) ,  stdout,printf  ( ),f close ( ) 

#include  "y_2d_ops . h"// . Rmatrix2D( ) , Rotate2D( ) 

int  main(){//<============================USE  ROTATIONS  TO  MODEL  CIRCULAR  ORBITS 


using  namespace  y2D0ps; 

FILE  *f=f reopen ("orbit .txt", nw",stdout);// . redirect  output  to  a  file 

double  s[2]={0,0},p[2]={l,0},m[2]={1.2J0};// . sun,  planet,  &  moon  positions 

double  Rp[4];/*<-*/RMatrix2D(Rp, .00436);// . rotation  matrix  for  planet 

double  Rm[4];/*<-*/RMatrix2D(Rm,  .279);// . rotation  matrix  for  moon 

printf("#  planet  |  moon\n#  x  ,  y  |  x  ,  y\n");//. column  header 

for(int  i=0;i<1444;++i){ 

printf  ("%6. 3f ,%6. 3f  ,%6. 3f  ,%6. 3f \n",p[0],p[l]  ,m[0],m[l] )’,//. .  record  positions 
Rotate2D(p,s,Rp),Rotate2D(m,s,Rp);//. .. .rotate  planet-moon  system  around  sun 

Rotate2D(m, p,  Rm) ; }// . rotate  moon  around  planet 

fclose(f ); 

}// - YAGENAUT@GMAIL.COM - LAST~UPDATED~02MAY2013 - 


orbit.txt 


#  planet  |  moon 

#  x  ,  y  |  x  ,  y 
1.000,  0.000,  1.200,  0.000 
1.000,  0.004,  1.192,  0.060 
1.000,  0.009,  1.169,  0.116 


EXAMPLE  2 


#include  <cstdio>// . FILE, freopen( ), stdout, printf ( ) , f close( ) 

#include  "y_2d_ops . h"// . <cmath>, IParameters2D( ) , Intersect2D( ) 

int  main(){//<======================USE  LINE-LINE  INTERSECTIONS  TO  MAKE  A  SPIRAL 

using  namespace  y2D0ps; 

FILE  *g=f reopen ("spiral .txt", "w", stdout);// . redirect  output  to  a  file 

double  X[250] [2]={  {-1,1}  ,  {-1,-1}  ,  {1,-1}  ,  {1,1}  }j 
for(int  i=3;i<249;++i){ 

double  a=(1.585)*(i-l)- .0285;// . "a"  increments  by  about  90  degrees 

double  A[4]={X[i] [0],X[i] [l],X[i] [0]+cos(a),X[i] [l]+sin(a)}; // . line  #1 

double  B[4]={X[i-3]  [0],X[i-3][l],X[i-2][0],X[i-2][l]};// . line  #2 

double  t[2] ;/*<-*/IParameters2D(t,A, B);// . calculate  intersect  parameters 

Intersect2D(X[i+l] , t, A) ; }// . use  parameters  to  find  X[i]  intersections 

printf("#  x  ,  y\n");// . print  column  header 

for (int  i=0;i<250;++i) printf ( "%6. 3f ,%6.3f\n",X[i] [0],X[i][l]);//. . print  points 
fclose(g) ; 

}// - YAGENAUT@GMAIL.COM - LAST~UPDATED~02MAY2013 - 


spiral.txt 


#  X  ,  y 
-1.000,  1.000 
-1.000,-1.000 
1.000,-1.000 
1.000,  1.000 
-1.000,  1.000 
-0.972,-1.000 
1.000,-0.944 
0.917,  1.000 


FIGURE  2 

image  created  from  spiral.txt 


. sin(),cos(),abs() 

TWO-DIMENSIONAL  OPERATIONS  (CARTESIAN  COORDINATES) 
;)@@@@@@@@@@@@@@@@@@@@@@@@@ 
PERFORMS  A  2D  TRANSLATION 
COORDINATES  TO  TRANSLATE  (MODIFIED  BY  THIS  FUNCTION) 
. DISPLACEMENT  VECTOR 


NO.  OF 

COPIES  ORGANIZATION 


1  DEFENSE  TECHNICAL 
(PDF)  INFORMATION  CTR 
DTIC  OCA 

1  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
IMAL  HRA 

1  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
RDRL  CIO  LL 

1  GOVT  PRINTG  OFC 
(PDF)  AMALHOTRA 

24  DIR  USARL 
(5  HC,  RDRL  WM 
19  PDF)  P  BAKER 
RDRL  WML 
M  ZOLTOSKI 
RDRL  WML  A 
M  ARTHUR 
B  BREECH 
P  BUTLER 
B  FLANDERS 
W  OBERLE 
C  PATTERSON 
R  PEARSON 
L  STROHM 
A  THOMPSON 
R  YAGER  (5  HC) 

RDRL  WML  B 
N  TRIVEDI 
RDRL  WML  C 
S  AUBERT 
RDRL  WML  D 
R  BEYER 
RDRL  WML  E 
P  WEINACHT 
RDRL  WML  F 
D  LYON 
RDRL  WML  G 
J  SOUTH 
RDRL  WML  H 
J  NEWILL 
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